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I. INTRODUCTION 



The radiative correction of the relative order a(Za) 3 to the energies of positronium 
and the radiative-recoil correction of the relative order a(Za) 3 ^ to the hydrogen energies 
should be taken into account when one compares the recent experimental results |l],|2| 
with the QED predictions. An attempt to calculate such the corrections, induced by 
the vacuum polarization, was undertaken in Ref. || for hydrogen. In Ref. the same 
correction was calculated for a two-body system of particles with an arbitrary mass ratio 
and the result was applied to hydrogen and positronium. 

In the present note the calculation of the vacuum polarization correction is reex- 
amined with the purpose to clarify the issue of the infrared regularization. Such a 
regularization is necessary for a control over the linear infrared divergence arising in the 
course of the calculation and through which the leading order correction reveals itself. 

In order to separate the contribution of the leading order, which comes from the 
atomic scale, from the next-to-leading one, saturated by the relativistic scale, we intro- 
duce the auxiliary parameter A', mZa C A' < m. Subtracting the Yukawa potential 
from the Coulomb one, 

Za Za Za w_ , , 
► + — e" Ar , 1 

ry ry ry 

we find the contribution to the energy due to the range of r's satisfying ^7 < r ~ On 
the other hand, the subtracted contribution is saturated by the short range, ^7 > r ~ — . 
In the sum of two contributions, the dependence on the parameter A' is cancelled away. 

Below we use a standard scheme of calculations, taking the vacuum polarization into 
account by means of the substitution 



(2) 



for a propagator of the photon, which polarizes vacuum. An energy correction which 
thus becomes a function of A is then integrated over A with a density of intermediate 
states p(A), which equals 



2a\ z + 2m 2 I 4m 2 
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for a free particle-antiparticle pair. 

II. ORDER ma(Z«) 4 CONTRIBUTION 

If bounded particles are non-relativistic, vacuum polarization manifests itself primar- 
ily due to the Coulomb field. Since the atomic momentum is much less than a mass A 
of the particle-antiparticle pair, this momentum can be neglected in the right-hand side 
of (0), so that an effective potential induced by the Coulomb vacuum polarization, at 
fixed A turns out to be 

Vx(r) = -^5(r), (4) 

while the corresponding (lowest-order) contribution to an atom's energy is 

47rZo#(0)| 2 

£lo(A) = ^ • (5) 

After the integration with respect to A with the weight function (|3|) it turns into 

AElo = jj^pi 

Now let us find the correction to this result induced by the modification of the 
Coulomb potential at short ranges (recall that A' ^> mZa) , 

Za Za Za y r 
-> + e~ Xr . (7) 
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To this end we calculate the average value of the operator (f|) over the state, whose wave 
function is perturbed by the short-range correction to the Coulomb potential (|7|): 

Zae- X< 



V x (r')g(r',r\E) 



ip(r) ) + c.c. 
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Here Q{r',r\E) is the reduced Green's function of the Schrodinger equation in the 
Coulomb field, ip{r) is the solution of this equation. To the lowest order in Za, we can 
disregard the Coulomb interaction in Q and the atomic momentum as compared with 
A': 



0(0,40)-- 
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^(r)-V(O), 
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(/i = mMj {m + M) is the reduced mass) and obtain 

An(Za) 2 



(10) 



II. ORDER ma{Zaf CONTRIBUTION 



In the next-to-leading order, let us consider the double photon exchange accounting 
for vacuum polarization by one of these photons. According to the Feynman's rules, 
corresponding contribution to the energy at fixed A is 

c ( \ \i\ 0/7 ^2 ,2 f dk o f d 3 k 4tt 4vt / (l + 7o)m + 7 /i;o-7fc \ 
£(X, A ) = -2(Za) * ]—J ^ 

(r /t k 2_Q2 T v -(ji^U,M^-M)J. (11) 

This expression can be represented graphically as two fermion lines connected by two 
photon ones. The first and the second terms in the angle brackets correspond to the 
graphs with uncrossed and crossed photon propagators, respectively. In (|TTD, k 2 = k 2 , 

K = ^k 2 - A 2 , K' = ^k 2 - A' 2 ; 
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lu = \jkl + 2mk , Q± = yJk%±2Mk ; 



Ifi (r^) are the Dirac matices for the light (heavy) particle. The parameter A will be 
used below as a mass of the virtual pair, in accord with (D), while the parameter A' is 
introduced to regularize the otherwise infrared divergent integral in (pT|). Two possible 
ways to insert A and A' into the photon propagators are accounted for in (|TT|) by the 
overall factor 2. Choosing A, A' ^> mZa, we can neglect atomic momenta so that taking 
the average over a bound state reduces to that over the Pauli spinors which is denoted 
by angle brackets, together with the multiplication by tp 2 = |^(0)| 2 . The spinor averages 
are trivial: 

(7m(i + 7bh*> (r^)(i + r )r^)) = 4, (12) 
(7m7o7^) + r )r„^)) = (7^(1 + 70)7^) (r^Tor^) = 2, (13) 
(7m7o7.) (r MM r r^)) = 4 T 20V, (14) 
(7^77,) (r^rr^) = -6 ± a** 1 . (15) 



Before proceeding further, let us consider the analytic properties of the integrand in 
(PI) as a function of ko- The photon propagators have the poles at the points ±\J k 2 + A 2 
and ±y k 2 + A' 2 . After the integration over k these poles turn into the cuts that go from 
— 00 to —A (—A') and from A (A') to oo. Similarly, the light fermion propagator gives 
rise to the cuts (— oo, — 2m] and [0, oo). The heavy fermion propagator in the second 
term of (|i~Tf) , which corresponds to the graph with crossed photon lines, gives rise to the 
cuts (— oo, — 2M\ and [0, oo). Finally, the heavy fermion propagator in the first term, 
corresponding to the graph with uncrossed photon lines, produces the cuts (— oo, 0] and 
[2M,oo). 

It is convenient to extract the heavy fermion propagator from the first term of (|iT|), 

1 1 _ 1 / 1 1 

k 2 -u 2 k 2 - Q 2 _ ~ 2(M + m)k [k 2 -u 2 ~ k 2 -Q 2 _ 
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and then to change the sign of the integration variable k Q — > — k in all terms containing 
l/(/c 2 — In this way the integral over k is naturally splitted into two parts. The 
former one, 

_ (Za) 2 ip 2 f dk n f d 3 k 4vr 4tt 



(2Mm 
\~kV 



M 2 - m 2 Jc- 2ni J (2tt) 3 k 2 - K! 1 k 2 - JC 



(17) 



+ 2m + 



m2k 2 -k 2 aa' 3k 2 - 2k 2 ^ 



M kn 



+ 



2M 



k 



k 2 -n 2 



(M <-> m) 



where Q = Q + , is taken over the contour C_, wrapping the left cut. The latter part is a 
residue at the k = pole, which appears in terms containing kQ 1 {k 2 — VL 2 _ ) _1 after the 
change k -^ — k : 



S pole (X, A ) - -2 M + m A + y ^^r" + 1 - — 



18) 



Let us begin with the latter contribution. In (|IBD , the first term in the brackets 
is singular when one of the lambdas approaches zero. This term is just the regulator 
contribution subtracted from the leading order correction. In fact, neglecting A' as 
compared with A in the sum A + A' above, we see that the infrared-singular term in 
(|18|) is compensated by the effect of the Coulomb potential modification (|H]). Hence, 
only two last terms in the right-hand side of Eq. ([l8|) comprise the genuine order (Za) 5 
contribution of the k = pole, so that we can safely set A' = in those terms: 



p v ; M + m A I 3 
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Likewise, we can set A' = in the left-cut contribution fllTf ) as far as this procedure 
does not spoil the infrared convergence of the integral. The integration over k gives for 
£c«t(A) = £ cu t(\, 0): 

4(Za) 2 ib 2 f „ f (2m m\ 1 
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Finally, integrating with respect to k and adding up the pole contribution we 
obtain: 

2tt era' 
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(21) 



Here 



^(x) = 0(1 -a) 



y/1 — .r 

X 



cos 1 x — 9(x — 1)^-^- cosh 1 x. 

x 



(22) 



From (|2"ID, the known results for the vacuum polarization contribution to the hyper- 
fine splitting in muonium |5],|6| and positronium |7[] can be obtained: 

+ /r7 f3 m/ l2 M 8, M 28 tt 2 \ ^/m 2 \] , , 



A<f = ^(^)-Bf. (24) 



Here 



& = (25) 
3 Mm V ; 



is the Fermi splitting with the anomalous magnetic moments omitted. 

For the spin-independent part of the correction to the energy levels of hydrogen and 
positronium, we have 

A(Za) 2 i) 2 [™ . , [ 7T 1 



Mm Ji HK ' \ 2x 3 x 2 x 3 

+ ^ + In 2x - 2x 2 In 2x + 2xVx 2 - 1 cosh -1 x\ (26) 

_ «(^«) 2 ^ /"47tt 2 70 /m\\ _ 
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and 



9 r— i x cosh 1 x 1 . . 

-2 In 2x + 2x 2 In 2x - 82VX 2 - 1 cosh" 1 x - 2 ^ (28) 



a(Za) 2 i) 2 49ir 2 _ 40 
^ ^2 I ^88 27 



(29) 



where 



a 2x 2 



P( x ) = o 4 — ^x 2 - 1. 

These results differ from those obtained in Refs. for hydrogen and in Ref. M for 
positronium. The error made in both works has the same origin - inaccurate treatment 
of the infrared divergence. In fact, the authors of Refs. ||f| do not introduce in (0) the 
parameter A', which, as we have seen above, regularizes the infrared divergence. Instead, 
they subtract from the integrand in (|TT|) its asymptotic value at small fc's. Giving the 
finite results of Refs. 0,|J, this last procedure cannot be correct, since those finite results 
arise as a difference between two divergent integrals. In contrast, the regularization 
procedure used in the present work deals with the well-defined finite expressions. 



Numerically, the correction (27) constitutes 0.647 kHz for the ground-state energy 
of hydrogen, while ( p9|) equals 46.7 kHz for the ground-state energy of positronium. In 
the case of hydrogen, the correction exceeds the uncertainty of the recent measurement 
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